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Abstract 

We simulate the central reactions of ^0^6+20]^^^ ^^Ar+^^Sc, ^^Ni+^^Ni, 
86Kr+93Nb, i29xe+ii8Sn, ^^Kr+^^'^Au and i^'^Au+^^'^Au at different incident en- 
ergies for different equations of state (EOS), binary cross sections and different 
widths of Gaussians. A rise and fall behaviour of the multiplicity of intermediate 
mass fragments (IMFs) is observed. The system size dependence of peak center- 
of-mass energy E™^^ and peak IMF multiplicity (N/a/f)™""^ is also studied, where 
it is observed that E™^^ follows a linear behaviour and (N/jv/f)™''^^ shows a power 
law dependence. A comparison between two clusterization methods, the minimum 
spanning tree and the minimum spanning tree method with binding energy check 



(MSTB) is also made. We find that MSTB method reduces the (N/mf)""''' espe- 
cially in heavy systems. The power law dependence is also observed for fragments 
of different sizes at E™^^ and power law parameter r is found to be close to unity 
in all cases except A"^"^. 

1 Introduction 

The behaviour of hot and dense nuclear matter at the extreme conditions of temperature 
and density is a question of keen interest. It can be studied with the help of heavy-ion 
reactions at intermediate energies. At high excitation energies, the colliding nuclei may 
break into several small and intermediate size fragments followed by a large number of 
nucleons [1-3]. A large number of experimental attempts had been carried out ranging 
from the evaporation of particles to the total disassembly of the dense matter. Besides 
these two extremes, one can also have a situation where excited matter breaks into several 
fragments. In the last few decades, several experimental groups have carried out a com- 
plete study of fragment formation with in detectors [4-10]. It is quite obvious from these 
studies that the fragments formed in heavy-ion collisions depend crucially on the bom- 
barding energy and impact parameter of the reaction [1-5] . Therefore, these experimental 
studies of fragmentation offer a unique opportunity to explore the mechanism behind the 
formation of the fragments. Moreover, one can also pin down the role of dynamics in 
fragment formation and their time scale. 

Recently, there has been increasing interest in the effects of reaction dynamics on 
the production of IMFs and light charged particles (LCPs, Z—1 or 2). Sisan et al. [6] 
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studied the emission of IMFs from central collisions of nearly symmetric systems using 
47r- Array set up where they found that the multiplicity of IMFs shows a rise and fall with 
increase in the beam energy. They observed that E^^^ (energy at which the maximum 
production of IMFs occurs) increases hnearly with the system mass whereas a power law 
(oc A"^) dependence was reported for peak multiplicity of IMFs with power factor r = 
0.7. Peaslee et al. [7], on the other hand, studied asymmetric system ^^Kr+^^^Au in the 
incident energy range from 35 to 400 MeV/nucleon and obtained an energy dependence 
of multifragmentation. Their findings revealed that fragment production increases up to 
100 MeV/nucleon and then decreases with increase in incident energy. De Souza et al. [8] 
studied the central collisions of ^^Ar+^^^Au from 35 to 120 MeV/nucleon and observed 
that IMF multiplicity shows a steady increase with increase in the incident energy. The 
IMF multiplicity decreases, however, when one moves from central to peripheral collisions. 
On the other hand, Tsang et al. [5], in their investigation of ^^^Au+^^^Au collisions at 
E/A = 100, 250, and 400 MeV, found the occurrence of peak multiplicity at lower energies 
for central coUisions whereas it is shifted to higher energies for peripheral collisions. Stone 
et al. [9] used a nearly symmetric system of ®^Kr+^^Nb from 35 to 95 MeV/nucleon to 
obtain IMF multiplicity distribution as a function of beam energy by selecting central 
events. Ogilvie et al. [10] also studied the multifragment decays of Au projectiles after 
collisions with C, Al, and Cu targets at the bombarding energy of 600 MeV/nucleon 
using ALADIN forward spectrometer at GSI, Darmstadt, with the beam accelerated by 
SIS synchrotron. They found that with increasing the violence of collision, the mean 
multiplicity of IMFs originating from projectile first increases to a maximum and then 
decreases again. 
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As mentioned earlier, Sisan et al. [6] reported that the peak multiphcity of IMFs as 
well as peak center-of-mass energy scale with the size of the system. In a recent commu- 
nication, Vermani and Puri [11] succeeded partially in explaining the above mentioned 
behaviour by using the quantum molecular dynamics (QMD) approach. We here plan 
to extend the above study by incorporating various model ingredients such as equation 
of state, nucleon-nucleon (nn) cross section, and Gaussian width. The role of different 
clusterization algorithms shall also be explored. We shall attempt to find out whether 
these ingredients have sizable effects. 

2 The Formalism 

2.1 Quantum Molecular dynamics (QMD) model 

We describe the time evolution of a heavy-ion reaction within the framework of Quan- 
tum Molecular Dynamics (QMD) model [1-3, 12] which is based on a molecular dynamics 
picture. The explicit two- and three-body interactions lead to the preservation of fluctua- 
tions and correlations that are important for N-body phenomena like multifragmentation. 
In QMD model each nucleon is represented by a Gaussian distribution whose centroid 
propagates with the classical equations of motion: 



dvi dH 
dt dpi ' 
dpi dH 



(1) 



(2) 



dt dvi ' 



where the Hamiltonian is given by 




(3) 
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with 




+ v- 



MDI 



(4) 



V^"'^ is the Skyrme force whereas V'^"^\ yYuk yMDi jggj^g^ respectively, the Coulomb, 
Yukawa and momentum dependent potentials. Yukawa term separates surface which 
also play role in low energy process hke fusion and cluster radioactivity [13,14]. The 
momentum-dependent part of the interaction acts strongly in the cases where the system 
is mildly excited [15, 16]. In this case, the MDI is reported to generate a lot more frag- 
ments compared to the static equation of state. For a detailed discussion of the different 
equations of state and MDI, the reader is referred to Refs. [4,15,16]. The relativistic 
effect does not play role in low incident energy of present interest. 

The phase space of the nucleons is stored at several time steps. The QMD model does 
not give any information about the fragments observed at the final stage of the reaction. 
In order to construct fragments from the present phase-space one needs the clusterization 
algorithms. We shall concentrate here on the MST and MSTB methods only. 

2.2 Different clusterization methods 
2.2.1 Minimum spanning tree (MST) method 

The widely used clusterization algorithm is the Minimum Spanning Tree (MST) method 
[12]. In MST method, two nucleons are allowed to share the same fragment if their 
centroids are closer than a distance rmm, 



ri - rjl < r, 



(5) 
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where ri and rj are the spatial positions of both nucleons. The value of r^m can vary 
between 2-4 fm. This method cannot address the question of time scale. This method 
gives a big fragment at high density which splits into several light and medium mass 
fragments after several hundred fm/c. This procedure gives same fragment pattern for 
times later than 200 fm/c, but cannot be used for earher times. 

2.2.2 Minimum spanning tree method with binding energy check (MSTB) 

This is an improved version of normal MST method. Firstly, the simulated phase-space 
is analyzed with MST method and pre-clusters are sorted out. Each of the pre-clusters is 
then subjected to binding energy check [11]: 

We take Ehind — -4.0 MeV if > 3 and -Ebind = 0.0 otherwise. Here is the number 
of nucleons in a fragment and is center-of-mass momentum of the fragment. This is 
known as Minimum Spanning Tree method with Binding energy check (MSTB) [11]. The 
fragments formed with the MSTB method are reliable and stable at early stages of the 
reactions. 

3 Results and Discussion 

We have simulated the central reactions of ^ONe+^ONe (E;„6 = 10-55 AMeV), ^^Ar+^^Sc 
(Eiab = 35-125 AMeV), ^^Ni+^^Ni (E,„6 = 35-105 AMeV), ^^Kr+^^Nb (Eiab = 35-95 
AMeV), i29xe+ii8Sn {Eiab= 45-140 AMeV), s^Kr+i^^Au {Eiab= 35-400 AMeV) and 
^^^Au-F^'^^Au (Eiab = 70-130 AMeV). The energies are guided by experiments [5-7]. For 

6 



the present study, we use hard (labeled as Hard), soft (Soft), Hard with MDI (HMD) and 
Soft with MDI (SMD) equation of state. We also use standard energy-dependent Cugnon 
cross section [16] and constant isotropic cross section of 55 mb strength in addition 

to two different widths of Gaussian L = 1.08 and 2.16 fm^ (L*"'"'^'^). The superscript to the 
labels represent cross section. The phase-space is clusterized using clusterization methods 
described previously. The reactions are followed till 200 fm/c but the conclusions do not 
change when the reaction is over employing the validity of both algorithms. 

In fig. 1, we display the time evolution of IMFs for the reaction ^^Kr-|-^^Nb at incident 
energy of 75 AMeV employing MST method. In fig. 1(a), we display the model calcula- 
tions using Hard*^"^ (solid line) and Soft™^ (dashed hue). Prom fig. 1(a), we see that the 
number of IMFs are larger in case of Soft as compared to Hard. This is because of the 
fact that soft matter will be easily compressed as a result of which density achieved will 
be more which in turn will lead to the large number of IMFs as compared to that in case 
of Hard. It is worth mentioning here that the effect could be opposite at higher energies. 
Since at higher energies the IMFs may further break into LCPs and free nucleons. In fig. 
1(b), we display the results for Hard'^'"^ and Hard^^ (dotted line). As evident from the fig. 
1(b), the number of IMFs are nearly same for both type of cross sections. This may be 
due to the fact that for the central collisions, since the excitation energy is already high 
therefore, different cross sections have a negligible role to play. In fig. 1(c), we display the 
results for Hard along with two different widths of Gaussian i.e. L and L''^°"'^ (dash-dotted 
line). We find that the width of Gaussian has a considerable impact on fragmentation. 
As we change the Gaussian width (L) from 4.33 to 8.66 fm^, the multiphcity of IMFs is 
reduced by ~ 30%. Interestingly, the kaon yield also get reduced by the same amount 
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Figure 1: The time evolution of IMFs (5<A <44) for the reaction of ^^Kr+^^Nb at incident 
energy of 75 AMeV for different model ingredients . 

[17]. Due to its large interaction range, an extended wave packet (i.e. L*''°"'^) will connect 
a large number of nucleons in a fragment, as a result it will generate heavier fragments as 
compared to one obtains with smaller width. It is worth mentioning here that the width 
of Gaussian has a considerable effect on the collective flow [17, 18] as well as pion produc- 
tion also [17, 19]. In fig. 1(d), we display the results using Hard and HMD (dash-dot-dot 
line). Again the number of IMFs are nearly same for both EOS. This is expected since 
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the effect of MDI will be small at these energies. However, the scenario is completely 
different at high energies. Since at high energies, due to the repulsive nature of MDI, 
there is large destruction of initial correlations and the additional momentum dependence 
further destroys the correlations reducing further the multiplicity of IMFs. This leads to 
the emission of lots of nucleons and LCPs [20] . 

In fig. 2, we display the average multiplicity of IMFs, (Nimf), as a function of incident 
energy in the center-of-mass frame (Ecm.) for ^^Ni+^^Ni reaction employing MST (open 
symbols) and MSTB (sohd symbols) methods. Figs. 2(a) and 2(b) are for Soff^"^ and 
Hard'^"^, respectively. Lines represent the quadratic fit to the model calculations. In both 
cases, the number of IMFs first increases with incident energy, attains a maxima and then 
decreases in agreement with the previous studies [5-7, 10, 11]. Clearly, (N/mf) is more for 
MST method as compared to MSTB method. Since in case of MSTB method along with 
spatial correlations, an additional check for binding energy is also used, therefore it filters 
out the loosely bound fragments which will decay later. Hence, the fragments obtained 
with MSTB method are properly bound. A similar trend is obtained for all other reactions 
as well as different model ingredients used in the present study but is less pronounced 
in lighter systems like ^°Ne+^°Ne, ^°Ar+^^Sc as compared to heavier systems. However 
for Gaussian width L'"'°"'^, the (N/mf) is nearly zero in this incident energy range using 
MSTB method (not shown here). This is due to the fact that an extended wave packet 
(i.e. L'""""'') connects a large number of nucleons in a fragment, as a result it generates 
heavier fragments and the additional binding energy check further excludes the unbound 
fragments. 

In fig. 3, we display the peak center-of-mass energy E^^^ as a function of combined 
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Figure 2: The mean IMF multiplicity, (N/jv/f), as a function of incident energy in center-of- 
mass frame, Ec.m., for the reaction of ^^Ni+^*Ni. Solid (dashed) curves show the quadratic 
fits to the model calculations for MSTB (MST) method to estimate the peak center-of- 
mass energy at which the maximal IMF emission occurs. 

mass of the system employing MST method. Lines represent linear fitting proportional 
to mA. We find that the mass dependence of E^^"^ is insensitive to different EOS (fig. 
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Figure 3: The E^^^ as a function of composite mass of the system (A). The different lines 
represent the linear fits. Comparison of model calculations is made with experimental 
data [5-7] (solid stars). The percolation calculations [6] (solid triangles) are also shown 
in figure. 

3a), nn cross section (fig. 3b) as well as the width of Gaussian also (fig. 3c). It is slightly 
sensitive to MDl because for heavy systems EJT'^'^ is more as a result of which the effect 
of MDl becomes non- negligible. In fig. 3, the model calculations are also compared with 
experimental data [5-7] . It is clear from the fig. 3 that model calculations for E"^^ agree 
with experimental data [5-7]. This behavior is consistent for all the different choices of 



11 




90 180 270 360 450 90 180 270 360 450 
System mass (A) 



Figure 4: (^imf)"^"'^ as a function of composite mass of the system (A). The different 
hues represent the power law fits (oc A^). Comparison of model calculations is made with 
experimental data [5-7] (solid stars) . The percolation calculations [6] (solid triangles) are 
also shown in figure. 

model ingredients. 

In fig. 4, we display the peak multiplicity of IMFs {^imf)"^"'^ as a function of com- 
bined mass of the system employing MST method. The lines represent power law fitting 
proportional to A"^. The multiplicity of IMFs, in case of ^"^Ne+^^Ne and ^"^Ar+^^Sc, is 
obtained by excluding the largest and second largest fragment, respectively, to get the ac- 
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curate information about the system size dependence. (Nimf)"^"'^ sue obtained by making 
a quadratic fit to tfie model calculations for (N/mf) as a function of (Ec.m.). We find that 
the peak multiplicity is insensitive to cross section (fig. 4b) and MDI (fig. 4d) (for expla- 
nation see discussion of fig. 1). It is slightly sensitive to EOS (fig. 4a but highly sensitive 
to the Gaussian width (fig. 4c). On increasing the width of Gaussian, (N/mf)"*"^ reduces 
to a large extent. As discussed earlier, an extended wave packet (i.e. L^*""""^) will connect 
a large number of nucleons in a fragment, as a result it generates heavier fragments as 
compared to one obtains with smaller width. From fig. 3, we see that E^^^ shows linear 
dependence (oc mA) whereas (Nimf)'^"'^ (fig- 4) follows power law behaviour (oc A"^) with 
T nearly equal to unity. In fig. 4, the model calculations are also compared with experi- 
mental data [5-7]. It is clear from the fig. 4 that, as the system mass increases difference 
between model calculations and experimental results goes on increasing. This behavior 
is consistent for all the different choices of model ingredients. This may be due to the 
fact that the fragments obtained with MST method are not reliable because this method 
makes sense only when matter is diluted and well separated. This is true only in case of 
high beam energy and in central collisions. Therefore, we have to look for other methods 
of clusterization. As mentioned earlier, the fragments obtained with MSTB method are 
properly bound and reliable. So, as a next step, we check system size dependence of 
and (Nimf)"^""^ by using MSTB method for clusterization. 

In fig. 5, we display the E"^^ (left panels) and (N/mf)"'"'' (right panels) for Soft^"» 
(upper panels) and Hard'^"^ (bottom panels) as a function of combined mass of the system. 
Solid (open) symbols represent MSTB (MST) method. Prom left panels we find that 
^c!m. remains insensitive to the choice of clusterization method. The same is true for 

13 



35 
30 
25 

^ 20 
g 15 
g 10 



> 

CD 



5 

^ 30 
t . 25 
20 
15 
10 
5 



J 1 1 1 1 M 1 

Lsofr^ 


M 1 M 1 1 1 1 1 1 1 M 1 1 M IL 






E t 




r <°x' 
[mii|iii 


Data [j-/J i 
^ Percolation I 


= Hard'" 








- 

- ^./J 






III 1 1 " 



1 1 1 1 I I I 1 1 I I I I I I 1 1 1 1 I I I 1 1 I I I I I 
-■ ■ MSTB ° - 




20 
15 
10 
5 



90 180 270 360 450 90 180 270 360 450 

System mass (A) 



Figure 5: The E^^^ (left panels) and (N/mf)™'"^ (right panels) as a function of composite 
mass of the system (A) using Soft"^"^ (upper panels) and Hard"^"^ (lower panels) employing 
MSTB and MST methods. The different lines in left (right) panels represent the linear 
fits (power law fits). Comparison of model calculations is made with experimental data 
[5-7] (solid stars). 

(N/mf)™"""^ (right panels) but in low mass region. As the system mass increases, the 
(N/mf)™""^ becomes more and more sensitive to the method of clusterizaton. The MSTB 
method excludes the loosely bound fragments thus reducing the peak IMF multiplicity. 
The effect is uniform for both the EOS as well as for different cross section (not shown 
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Figure 6: The largest fragment and multiplicities of free-nucleons, LCPs, MMFs, HMFs, 
and IMFs as a function of composite mass of the colliding nuclei (A) for different reactions 
at their respective E™^^ (solid circles). Lines represent the power law fits (oc A^). 

here). 

In fig. 6, we display peak multiplicity (obtained by employing MSTB method) as a 
function of composite mass of the system for various fragments consisting of the largest 
fragment (A™"^) (fig. 6a), free-nucleons (1<A <1) (fig. 6b), light charged particles 
(LCPs) (2<A <4) (fig. 6c), medium mass fragments (MMFs) (5<A <9) (fig. 6d), heavy 
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mass fragments (HMFs) (10<A <44) (fig. 6e) and intermediate mass fragments (IMFs) 
(5<A <44) (fig. 6f) for Hard™^. Lines represent the power law fitting proportional to 
A"^. Interestingly, the peak multiplicities of different fragments follow a power law (oc A"^). 
Power law factor r is almost unity in all cases except A"*"^ for which there is no clear 
system size dependence. The system size dependence of various fragments has also been 
predicted by Dhawan and Puri [21]. Their calculations at the energy of vanishing flow 
(i.e., the energy at which the transverse flow vanishes) clearly suggested the existence of 
a power law system mass dependence for various fragment multiplicities. 



4 Summary 

We have simulated the central reactions of nearly symmetric, and asymmetric systems 
over the entire periodic table at different incident energies for the different equations 
of state (EOS), nn cross sections and different widths of Gaussians. We have observed 
that the multiplicity of intermediate mass fragments (IMFs) (3 <Z< 20) shows a rise 
and fall with increase in beam energy in the center-of-mass frame as already predicted 
experimentally/theoretically. We have also studied the system size dependence of peak 
center-of-mass energy E'^^ and peak IMF multiphcity (N/mf)"^"^- It has been observed 
that E™^^ increases linearly with system mass whereas a power law (oc A"^) dependence has 
been observed for (N/mf)"*"^ with r ~1.0. We have compared system size dependence 
of E"^^ and (N/mf)'""'' for MST and MSTB methods and found that MSTB method 
reduces the (N/mf)'""^ especially in heavy systems because in MSTB method due to 
binding energy check loosely bound fragments get excluded. The power law dependence 
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is also observed for fragments of different sizes at the energy for which the production of 
IMFs is maximum and power law parameter r is found to be close to unity in all cases 
except A"^"^. 
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